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INTRODUCTION 

In the past one of NASA's primary emphasis has been on identifying single 
and multiphase flow experiments which can produce new discoveries that are not 
possible except in a microgravity environment. While such experiments are 
obviously of great scientific interest, they do not necessarily provide NASA with 
the ability to use multiphase processes for power production and/or utilization in 
space. 

The purpose of the research presented in this paper is to demonstrate the 
ability of multidimensional two-fluid models for bubbly two-phase flow to 
accurately predict lateral phase distribution phenomena in microgravjty 
environments. If successful, this research should provide NASA with 
mechanistically-based analytical methods which can be used for multiphase 
space system design and evaluation, and should be the basis for future shuttle 
experiments for model verification. 

DISCUSSION 

During the last decade mechanistically-based multidimensional two-fluid 

models have been developed and successfully applied to the prediction of bubbly 

two-phase flows. It appears that these models should also work in microgravity 

envir onmen ts, however this s till needs to be verified. 
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To this end experiments are being conducted within Rensselaer's Center for 
Multiphase Research (CMR). In particular, microgravity conditions have been 
simulated using neutral buoyant polystyrene spheres, and future experiments 
will be performed in which neutral buoyant oil droplets are immersed in flowing 
water. 

The purpose of this paper is to present the progress to date in the analytical 
modeling of dispersed flows and to present the solid/fluid data which has been 
taken to support these modeling efforts. 

We will begin by summarizing the multidimensional two-fluid model and then 
will present the solid/fluid data which has been acquired in this program. 

ANALYSIS 

The analysis of multidimensional two-phase flow can be done using two-fluid 
models and associated computational fluid dynamic (CFD) numerical evaluation 
algorithms (eg, PHOENICS or FLOW3D). For example, the evaluation of 
adiabatic bubbly two-phase flows are governed by the two-fluid conservation 
equations for mass and momentum. These balance equations can be derived 
using ensemble averaging techniques [Lahey & Drew, 1992]. 

T WQ - P H AS E F L OW BAL A NC E EQUATIO N S 

The three dimensional balance equations for adiabatic two-phase flow are: 

Mass . Conser va tio n 
d[ a kPkl 

— ^+V.[a k p k v k ] = 0 (1) 

where a k is the volume fraction, p k is the density, and v^ is the average velocity of 
phase-k, respectively. 
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5[<XkPkXkl r Re 1 

+ V*[a k p k Ykv k ] = - V[a k p k ] + V* [a k (jk + ^ )J + a k Pk| + Mk 


(2) 


where M. k is the interfacial force density, p k is the static pressure, x^ is the 

Re 

viscous stress tensor, and Xk is the Reynolds stress on phase-k, respectively. 

The interfacial jump condition for momentum is: 

+ (3) 

where, m? is the interfacial momentum source density due to surface tension 
effects (eg, Marangoni forces). 


CLOSURE 

For dispersed spherical particles (eg, bubbles), we may use the interfacial 

transfers developed by Park [1992]. It is conventional to partition the interfacial 
force density (M^) into drag (d) and nondrag (nd) components: 


,, „(d) __(nd) 

M v = My + M v • 


We assume the following form of the interfacial drag law: 
M^ d) = - M^ d) = | P^ C D ly r lv r A'" 


(4) 

(5) 


where, y r = y v - y.£, and for spherical monodispersed bubbles, A"' = 6a y /D b is the 

interfacial area density of the interface between the continuous phase and the 
dispersed bubbles. Assuming the validity of inviscid flow theory for the 
continuous phase the non-drag interfacial force density for the continuous phase 
is, using cell-model averaging techniques [Park, 1992], [Lopez de Bertodano, 1992]: 




= - P^.Vcv + o^C^aym + OvP^C^-Vy,. 
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+ a v P i C 2 [x r *(V 2 r') + (V*Y r )v r J + bgP^r-YrVOv 

+ a s p^ r *Vot v Y r + OyP^C^Xr X V X Y v + OyPfi^r xVx 2< ( 6 ) 

Also, for a spherical dispersed phase the interfacial force required to have the 
bubbles maintain a spherical shape is [Park, 1992]: 

BLi = V«[a y Og] = V^a v p^ (a s v r v r +b s y r *v r l)J (7) 

Also, for the interfacial averaged- pressure, we have [Park, 1992]: 


P/.-Pz + CpP/ 1 ^ 


(8) 


For bubbly two-phase flows the total Reynolds stress tensor for the continuous 
liquid phase is given by superposition as, 


Re Re Re 
Ji =j£(BI) + J*(SD 


(9) 


where, for bubbly two-phase flows the bubble-induced shear stress is given by cell- 
model averaging [Lopez de Bertodano, 1992] as: 

Re 


«^(BD = * W [ a ^r + Vr^rl] 


( 10 ) 


We note that^gj) is the shear- induced Reynolds stress which may come from a 


classical k-e model and an algebraic stress law [Rodi, 1984]. 

For the inertial coupling of all dispersed/continuous phase interactions we 
have [Park, 1992]: 


Cym = ^T. + = O » ^T. - ^rot. “ A > ^1 “ A > 


'rot ” 2 ’ '“'L - '-'rot “ 4 » '-'l “ 4 


r_— - _ _ 9 _ , £ __ 3 _ p_l 1 

U 2 = - 20 . a s “ a s “ " 20 ’ D s =D s “ 20 ’ U P “ 4 ’ a i ~ " 20 
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It is significant to note that there are no arbitrary constants in the two-fluid 
model, however data on the interfacial drag coefficient, C D , and the lateral lift 

coefficient, C L , indicate that they should be a function of Reynolds number. 

MODEL ASSESSMENT 

The two-fluid model given in Eqs. (l)-(ll) has been assessed against a wide 
variety of terrestrial air/water bubbly flow data. 

Figure- 1 shows a comparison of the two-fluid model with the air/water bubbly 
upfiow void fraction data of Serizawa [1974]. Good agreement can be seen. 
Figure-2 shows that the turbulence model being used does a good job of predicting 
the measured nonisotropic turbulent structure of Serizawa. Similarly, Figure-3 

shows that the model predicts Serizawa's two-phase Reynolds stress data and 
that, unlike fully developed single-phase flow (ie, j = 0), the two-phase Reynolds 

stress distribution is not linear in the radial direction (since the radial void 
distribution is not linear). 

Figure-4 is very exciting because it shows that the same multidimensional 
two-fluid model which was used to predict Serizawa's bubbly upflow data is able to 
predict the air/water bubbly downflow data of Wang et al [1987]. Notice the lateral 
phase distribution in Figure-4 is completely different from Figure-1. This is due 
to the fact that the lateral lift force changes sign for downflows since the effect of 
buoyancy causes the relative velocity to change sign. 

Figures-5&6 are perhaps the most impressive comparisons of all, since they 
show that the same multidimensional two-fluid model is also able to predict the 
lateral phase and velocity distribution in complex geometry conduits (ie, vertical 
air/water bubbly upflows in an isosceles triangular test section). 

These data comparisons clearly indicate that the multidimensional two-fluid 
model given in Eqs. (l)-(ll) is able to predict a wide range of adiabatic bubbly flow 
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data taken on earth. Moreover, since this model is based on first principles it 
should also be able to predict lateral phase distribution for bubbly flows in 
microgravity experiments. 

EXPERIMENTS 

In order to simulate microgravity bubbly flow data, a experiment has been 
performed using approximately 2 mm diameter expanded polystyrene spheres 
immersed in water. The specific gravity of these spheres was 1.03, thus they were 
essentially neutral buoyant. 

Figure-7 is a schematic of the test loop. Figure-8 shows a schematic of the 
novel ventri-type phase separation device which was developed and used to avoid 
damage to the dispersed particles. Figure-9 is a schematic of the horizontal test 
section and the DANTEC three-dimensional laser Doppler anemometer (LDA) 
system which was used. It can be noted that the fiber optic laser 
transmitter/receiver heads were submerged in water and the tubular test section 
was constructed of fluorinated ethylene propylene (FEP) which has about the 
same index of refraction as water, thus no corrections for laser beam refraction 
were necessary. 

Figure- 10 presents the measured particle volume fraction distribution. These 
data were corrected for beam interruptions, etc. using the method of Alajbegovid, 
et al [1994]. It should be noted when these data were integrated across the cross- 
section they agreed to within ±1% with the corresponding global volume fraction 
data, which were taken using quick closing valves (Figure-11). 

Figures-12 & 13 show the Reynolds stress data for the dispersed particles and 
the continuous liquid phase (water), respectively, and figures-14 & 15 show data 
for the axial velocity fluctuations of the particles and the water, respectively. It 
can be seen that both data sets yield similar results. 
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Figure- 16 gives data on the mean axial velocity of the particles and the liquid 
phase, respectively. It can be seen that these data are symmetric (as they should 
be for neutral buoyant particles) and that the relative velocity is negligible, as 
would be expected for bubbly flow in a microgravity environment. Hence, except 
for interfacial boundary condition differences between a bubble and a solid sphere, 
these data comprise an excellent basis for the assessment of the two-fluid model 
for use in microgravity environments. 

MODEL ASSESSMENT 

Figures- 17 show good agreement between the multidimensional two-fluid 
model presented in Eqs. (l)-(ll) and the solid/fluid data discussed herein. 

Interestingly, it appears that the two-fluid model predictions agree better with 
the uncorrected particle volume fraction data than the corrected data. It is not 
completely clear at this time why this occurs, however the data correction method 
which we used [Alajbegovid et al, 1994] implicitly assumes that the dispersed 
particles are opaque, while the actual particles were translucent, thus it is likely 
that the data correction applied was inappropriate. This issue will require 
further study and the results of this study will be reported subsequently. 

Nevertheless, it appears that the essential physics of lateral phase separation 
is captured by the two-fluid model, and thus it should be appropriate for 
microgravity environments. 

Since solid particles satisfy the no slip condition at their surface while vapor 
bubbles and liquid droplets do not, a series of experiments using neutral buoyant 
oil droplets immersed in water will also be performed and the two-fluid model will 
also be assessed against these data. An oil/water loop, shown schematically in 
Figure- 18, has been designed and built and system shakedown is under way. 
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SUMMARY & CONCLUSIONS 

It has been shown that a mechanistically-based, multidimensional two-fluid 
model has been developed, and that it is capable of predicting a wide range of 
terrestrial bubbly flow data. 

An experiment has been performed in which detailed multidimensional 
measurements have been made for turbulent solid/fluid flows which simulate 
bubbly flows in a microgravity environment (ie, the spherical particles used had 
essentially the same density as water and thus the buoyancy term was 
eliminated). 

Initial comparisons between the solid/fluid data and the multidimensional 
two-fluid model showed good agreement and imply the validity of this model for 
microgravity environments. 
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particles with density equal to1030 Kg/m . 


Figure - 10 
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